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Abstract 

We study at the microscopic level the dynamics of a one-dimensional grav- 
itationally interacting sticky gas. Initially, N identical particles of mass m 
with uncorrelated, randomly distributed velocities fill homogeneously a finite 
region of space. It is proved that at a characteristic time a single macroscopic 
mass is formed with certainty, surrounded by a dust of non extensive frag- 
ments. In the continuum limit this corresponds to a single shock creating a 
singular mass density. The statistics of the remaining fragments obeys the 
Poisson law at all times following the shock. Numerical simulations indicate 
that up to the moment of macroscopic aggregation the system remains inter- 
nally homogeneous. At the short time scale a rapid decrease in the kinetic 
energy is observed, accompanied by the formation of a number ~ \fN of 
aggregates with masses ~ m\J~N. 
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1. INTRODUCTION 



The model of gravitationally interacting sticky particles has been pro- 
posed by Zeldovich [1], and then extensively studied in connection with the 
problem of large scale structures in the universe (for a recent review see 
[2]). The effort has been concentrated on the study of solutions of the cor- 
responding self-consistent hydrodynamic equations in order to understand 
how small density fluctuations arround a homogeneous state could induce 
observed mass distributions. 

In this paper we analyze the dynamics of a one-dimensional system with 
an initial homogeneous mass distribution filling only a finite region of space. 
As the one-dimensional gravitational interaction is confining, the whole sys- 
tem will eventually form a single mass. We start from the microscopic dy- 
namics of N point particles with randomly chosen initial velocities. Then we 
determine the evolution of the statistical distribution of masses formed by 
merging at binary sticky collisions. It turns out that a macroscopic mass is 
formed after a characteristic finite time with probability one, surrounded by 
a cloud of non-extensive fragments. In the continuum limit, this corresponds 
to a solution of the system of mass and momentum conservation laws, show- 
ing a single shock at the characteristic time. We emphasize that energy is 
not conserved in this process. We thus find not only a simple realization of 
a particular solution belonging to the general class studied in [3] , but also a 
detailed description of the statistics of masses in the course of time. 

As mentioned above, the system is composed of point particles attract- 
ing each other with forces proportional to the product of their masses, and 
independent of the interparticle distance. The corresponding n-body Hamil- 
tonian has the form 

n ^2 n n 

i={ 1 i < j 

where 7 is the gravitational constant, and rrii,Xi,pi denote the mass, the 
position and the momentum of particle i, respectively. 

The Hamiltonian ([!]) with an appropriate value of n defines the dynamics 
of the system in the time intervals separating binary collisions. The collisions, 
which are supposed to be perfectly inelastic, are responsible for the mass 
aggregation. When two neighbouring particles i and j meet, they merge 
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instantaneously forming a single mass (rrii + rrij) which continues the motion 
acquiring at the moment of impact the momentum (pi + Pj)- The complete 
dynamics, involving the aggregation process, is thus subject to the mass and 
momentum conservation laws. On the other hand, the number of particles 
decreases monotonically in the course of time, each collision replacing two 
particles by one. 

In a previous work [4] we concentrated on the determination of the prob- 
ability P/v(t) of merging before time t of the initial iV-particle dust into a 
single body (we denote here by iV the number of particles at time t = 0). 
The rigorous results have been obtained by assuming a uniform equidistant 
configuration of the initial identical masses m, characterized by a constant 
mass density 

m 

P=- 2 
a 

where a is the distance between the nearest neighbours. The velocities of the 
particles were supposed to be uncorrelated at t — 0, distributed according to 
some probability density 0(f). In particular, the Gaussian law 

has been considered in a detailed way. 

It has been found that for a macroscopic amount of matter N — > oo the 
probability -Pjv(t) of such a complete merging vanished for times shorter that 
the characteristic time 

** = 4= ( 4 ) 



IP 

More precisely, expressing the probability Pjv(i) in terms of the relevant 
time variable 

r = It - - (5) 
pt 

we derived a general (independent of the form of distribution 0), remarkably 
simple result 

Pn(t = 0) = 1 (6) 
from which the above assertion followed. 
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The exact form of the monotonically increasing function Pjv(t) was then 
determined in the limit of a continuous initial mass distribution. One way of 
defining this limit is given below: 

N -> oo, m -> 0, a -> (7) 

Mtot = Am = const, p = — = const, A = const 

a 

However, it has been stressed in our concluding remarks that the really 
relevant question was that of the macroscopic aggregation, which is quite 
different from the complete merging into a single mass. Indeed, from the 
physical point of view the important event is certainly that of the forma- 
tion of a macroscopic mass, representing a finite fraction of the total mass 
M to t = Nm in the continuum limit (0). And this event, according to nu- 
merical simulations and preliminary analytical results, seemed to occur in 
the immediate vicinity of the characteristic time t*, which could not be pre- 
dicted from the knowledge of the probability P^(t). The main object of the 
present study is to clarify this point by analyzing the evolution of the mass 
distribution in the course of time. In other words, we determine here the 
statistics of the mass aggregation providing a proof that the mass density in 
the continuum limit (|7|) becomes singular after a finite time t*, and takes the 
form of the Dirac 5 centered on a position of a single macroscopic mass. 

The concentration of the mass density on a single point at time t* can be 
simply demonstrated in the special case of a static initial condition, when all 
the particles at t = are at rest: <f>(v) = S(v). Indeed, suppose that the initial 

positions of the particles are Xj(0) = ja, j = 1,2, , N. Then, at time 

t > 0, the distances between neighbouring particles shrink to (a — rrryt 2 ), so 

that all the particles merge simultaneously at the moment t = t* = ^Ja/^m. 
The initial mass density uniformly distributed within the interval < x < Na 
is given by 

p(x; 0) = 9(x)6(Na - x)- = 6(x)6 ( — - x ) p, 

a \ p J 

where 9(x) is the Heaviside unit step function. For times < t < t*, it 
acquires in the continuum limit (^) the value 



p(x;t) = 9 



t 2 

x - ^M tot — 







M tot t 2 

jM tot — - x 

p 2 



Pit 2 
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Equation (|8|) implies the expected result 



lim p(x; t) = M tot 8 ( x 



tot 



(9) 



2p 



It is shown in the present paper how the above result can be generalized to 
the case of a random (Gaussian) initial velocity distribution (^), reflecting 
the creation of a macroscopic body at t — t*. Moreover, we derive an ana- 
lytic formula for the distribution of the remaining nonextensive microscopic 
masses. Computer simulations supplement our analytic approach on two 
points. They permit to determine the scaling laws to the continuum limit 
and to predict the mass distribution before t*, a problem which remains open 
for the mathematical analysis. 

Before closing these introductory remarks let us recall that the system 
of mass and momentum conservation laws governing the dynamics of one- 
dimensional aggregation has been recently studied from the point of view of 
the existence of global weak solutions [3]. In their proof the authors have 
analyzed the dynamics in terms of the center of mass trajectories in much the 
same way as it had already been done in our previously published papers [5,4]. 



2. EVALUATING THE MASS DENSITY 

At the initial moment t = 0, N identical masses m start the motion with 
uncorrelated velocities, distributed according to the probability density 4>(v). 
Particle j begins to move from the point ja, j = 1,2,...,N. In order to 
study the evolution of the mass density one has to determine the probability 
density for finding at time t > a mass M at point X. The dynamics 
of the system implies that M results from aggregation of some cluster of 
neighbouring initial masses, and that X must coincide with the position of 
the center of mass of the cluster at time t. So, let us consider a n-particle 
cluster 



In accordance with the dynamics induced by the Hamiltonian (|l|) its 





(10) 
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center of mass X™ +1 follows the trajectory 



with velocity 



3 + 



(n+ 1) 



t n 
a H — v 



n 



i+s 



r jm{N 



n 



=1 



2 % 



V; +1 (t) = -J2v J+s + 1 m(N 



n 



2j)t 



(12) 



8=1 



if 



The n particles (0) merge into a single mass nm before time t if and only 



Xr +1 (t)>X^ r r +1 (t), r = l,2,...,n-l (13) 

The inequalities ( |I~3"D express the requirement that for any partition of the 
n-particle cluster fllQl) into subclusters (j + 1, j + 2, j + r) and (j + r + 1, j + 
r + 2, n) the centers of mass of the subclusters cross before time t leading 
to the total merging. In order to guarantee that a single mass M n = nm is 
actually observed at the point ( [TT]) we have still to rule out the possibility of 
disturbance which would cause collisions with surrounding masses, built up 
from the initial clusters 

(14) 



:i,2, 



and 



(n+l,n + 2,....,N) 



(15) 



The unperturbed aggregation of the n-particle cluster ( |i~0D occurs if and only 
if 



X 



j-s+l 



(t)<X? +1 (t), 



1,2,...,J 



X? +1 (t) < X^ +1 (t), 8 = l,2,...,N-j 



n 



(16) 
(17) 



We continue to use here the notation X™ +1 to denote the position of the center 
of mass of the n-particle cluster composed of particles (i + 1, i + 2...., i + n). 

The inequalities fll6| ) express the fact that the center of mass trajectories 
of the s-particle clusters 



l,j -s + 2,...,j), 



1,2, 



■J 



stay to the left of the trajectory X™ +1 (t) up to time t, and thus do not cross 
it. Similarly, the inequalities ( 17|) exclude crossing of the trajectory of the 
aggregating mass M n with the center of mass trajectories on which evolve 
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masses ( |I5D initially to the right of it. The necessary and sufficient character 
of the conditions flUf ) and ( |17D follows from the remark that the dynamics 
excludes more than one crossing between the particle trajectories. 

The probability density for finding at time t > a mass M at point X 
with velocity V can now be written in the form 



N N-n n-1 

<EE^- Xj+iWW - V? +1 {t)]5{M - nm) J[ 6[X] +l {t) - X^ +l {t)} 

n=l j=0 r=l 
j N-j-n 

xn^+i(*)--x-/-^i(*)] n e[x^ +1 (t) - x? +1 (t)\ > (is) 

s=l s=l 

where < ... > denotes the mean value with respect to the initial velocity 
distribution 

N 

(19) 

8=1 

Multiplying formula ([18]) by M, and integrating over all possible masses and 
velocities we arrive at the expression for the mass density 

N N-n n-1 

n=l j=0 r=l 
j N-j-n 

x < J] B[X - X^ s+1 (t)\ >< I] e\X s n+1 {t)-X}> (20) 

s=l s=l 

In writing equation (^) we took into account the absence of correlations 
in the velocity distribution (plj). As a result the mean value appearing in 
equation (p~8| ) factorized out into the product of three averages correspond- 
ing to disjoint groups of velocity variables (vi,...,Vj), (vj+i, Vj +n ) and 
(vj +n+ i, wn). The calculations greatly simplify in the case of the Gaussian 
form (|3]) which will be used in the sequel. In order to show the kind of prob- 
lems one has to deal with let us consider in equation ( [2"0"D the integration 
over the velocities (vj + %, Vj + %, Vj +n ) of the aggregating n-particle cluster. 
Using equation ([□]) we find 

n-l 

< s[x - x? +1 {t)] n e[jq +1 (t) - x^; +1 (t)} > 

r=l 



s 



dv 



3+1- 



dv j+n <f) X {v j+1 )...<f>x (v j+n ) 



X-(j+ V -^)a - - £ v j+s - im {N -n- 2j)- 



n 



n-X 



X 



n 9 \ z) v j+s - - £ + 



r(n — r) 



-mr 



n 



(21) 



By introducing new integration variables 
1 



u. 



A 



\s=X 



n 



^2v j+s \ , r= l,2,...,n-l 



1 n 



A ^ ^ 



(22) 



one eventually rewrites formula ([21]) in the form 

l 



At V 2tt 
Here 



cxp 



n 




2A 2 





Mtot + m 
2p 





\N -n 




— mr 


[ 2 J _ 


)'] 



P„(t) (23) 



P n (r) = v / 2vm y dux— J du n -x<t>{u x )(l>{u2 - u x )-(j){u n ^x ~ u n _ 2 )<t>{-Un-x) 

r(n — r)mr 



n-X 
r=X 



U r + 



2A 



(24) 



is the probability of merging of an isolated n-particle cluster into a single 
mass before time t, with 



0H = ^exp UL 



(25) 



It is exactly this quantity which was the main object of our previous study 
[4]. In particular the important relation (^) has been derived therein. 

In quite a similar way one can analyze the other two averages appearing in 
equation (|20|) . We give the results here only for t = t* (or r = 0), as our aim 
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is to find the mass density at this characteristic moment. A straightforward 
calculation yields then the formulae 



3 

E 

s=l 



< n - *;-*x(**)] > = b 3 [ w 



X 



M tot + m 



N-j-n 

< II 0[Xn + l(t*) - X] > = B N ^ 
s=l 



1 



X- 



2p 

M tot + m 
2p 



(26) 
(27) 



where 



Bj(u) = J dv,\... J duj(j){u\)4>{u2 — U\)...(j)(uj — Uj-i) EI + u s ) 

= %[u s > -su, s = (28) 

As in [4], u s is interpreted as a Brownian path at discrete times s — 1,2, ... , 
and the traditional notation for the Wiener measure has also been 

used here. Taking into account the relation @ we eventually find that the 
mass density at the gravitationally imposed finite time scale t* in the con- 
tinuum limit (0) equals 



Y N N—n j 

p(X;t*) = M tot lim — V V — */— exp 

^ v ' ; n^oo N p At* V 2vr 1 



n 



2(At*) ; 



X - 



2p 



x b j < h— 

J V At* 



/<)/, 



2p 



5 



jV-n-j 



1 

At* 



X 



Mi 



lot 



2p 



(29) 



For symmetry reasons one can expect the macroscopic aggregation to 
take place at the central point of the originally uniform system. And indeed, 
it turns out that the mass density fl29|) becomes singular at this point in 
the continuum limit ([/]). In order to prove it we shall use now the Sparre 
Andersen theorem (Section XII. 7 in [6]) which permits to determine the 
generating function^ 

oo 

p(z;u) = l + Y, z n B n (u) (30) 



n=l 



1 B n (u) here is the same as B n in (42) of [4] with the variable u playing the role of r. 
Note that there is a minus sign missing in the equalities (46) and (48) of [4]. 
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\n=l 



oo z n 



n 



One finds (see the discussion in section 5 in [4]) 

/ 00 z n 

p(z; u) = exp ( ^2 — E w [u n > —nu] 

It follows that 

p(z; u)p(z; —u) = exp ( ^ — {E w [u n > —nu] + E w [u n > nu\ 

1 



Vn=l 



n 



ex p E 



00 z n 



\n=l 
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1 - Z 



Moreover, it follows from (|32|) that (setting Bo(u) = 1 



A' 



Y J z n B n {u)Y,^B t {u) = E z N Y,Be(u)B N ^(-u) = - 
n=o e=o n=o e=o 



which implies 



N—n 

E B J (u)B N 

-n-j\ 

j=0 



Thus the equation (|29|) simplifies in a remarkable way to the form 



1 N 1 rn~ 

t .) = M« i lim-i;-^/-exp 



n 



2(Xt*) 



X- 



Mi 



tot 



2p 



(31) 



(32) 



(33) 



(34) 



The main conclusion from equation (|34D is that at the moment t* the mass 



density concentrates at the central point of the system X = M tot /2p. Indeed, 
if f(X) is a continuous fonction, 



f{X)p(X;t*)dX 



1 N r ( 
M tot Jim ^E/^/( 



\ 2p y/n , 



1 / F 2> 
exp 



^7T 



AW 



2p 
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leading to 



p(X;t*) = M tot 5\X-^j (35) 



In particular, the formula 



p (M*EL. t *) „ Mtot-^ J2^-\ffi N^oo (36) 
yy 2p ' ; N ^ At* V 2vr' v ; 



shows that the mass density at X = M to t/2p (il) tends to infinity as v N. 

Equation (^) is the main result of this Section. It tells us that even in a 
nonstatic case where the aggregating particles have initially a Gaussian ve- 
locity distribution the macroscopic mass is formed at the gravitational time 
scale t = t*. From this point of view the situation does not differ from the 
static problem discussed in Section 1. However, the relation (|5|) clearly shows 
that the probability of complete aggregation at t — t* is still zero. This puts 
forward the question of the distribution of those masses which have not joint 
the central macroscopic body at t — t*. The statistical distribution of this 
leftover dust, representing a nonextensive amount composed of microscopic 
fragments, is discussed in the next Section. 



3. STATISTICAL DISTRIBUTION OF MASSES 

We first write down the probability distribution 

p k (XxVinx, . . . , X k V k n k ; t) 

for finding, at time t, an ordered configuration of k aggregates at positions 
Xi, (Xi < X 2 < . . . < X k ) with velocities Vi and masses M n . = mrii. 

The aggregates originate from consecutive initial clusters made ofrii,...,n k 
initial masses m and are found at the center of mass of these clusters 

/ n . _|_ \\ y 

X n .(t) = (nx + n 2 + • • • + n^ x + -^-j a + -^t + 7mA n - (37) 

with velocities 

K i (t) = — + im/\ n i (38) 
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The relations ( ftTD and (^) generalize ([□]) and (|P2| ) to the splitting of the 
iV equidistant initial particles into k consecutive clusters (here, i = 1, . . . , k 
does not index lattice sites, but the clusters themselves). In ( |3"TD and 



V ni is the sum of the initial velocities of the particles belonging to the i th 



cluster and 



mAiii = m{n i+ i H h n k - n x n;_i) (39) 



is the difference of the masses to the right and to the left of this cluster 
responsible for the force acting on it. 

The distribution /i^ is obtained by averaging the set of kinematical con- 
straints required for the realization of the desired event 

^(X 1 V 1 n 1 ,...,X k V k n k ;t) = (40) 

(l[e(x i+1 - x,) j[[6(Xt - x n xt))5{v t - v ni (t))} n ®nM) 

\i=l i=l i=l / 

We always have Y^=i n i — N, and the quantity Q ni (t) represents the 
constraint (to be elaborated below) needed to ensure that the i th aggregate 
forms before time t. Both X n .(t), V ni (t) and Q ni {t) are known expressions of 
the initial velocities and (...) denotes as before the mean value with respect 
to the distribution (0), so ^ can be calculated in principle. Notice that 
we did not include further constraints saying that the particles belonging to 
the adjacent clusters i and i + 1 do not perturb each other before t. This is 
not needed since the very fact that all initial particles in the (i + l) th cluster 
are found right of those in the i th cluster at time t already implies that 
no particles of the two groups have collisioned before t (otherwise a double 
crossing between trajectories would have occured, which is not possible in 
our dynamics). 

To simplify the discussion, we shall merely be interested in the mass 
distribution by integrating out positions and velocities 

H^(n 1 ,n 2 ,...,n k ;t) 
= JdX 1 ...dX k JdV 1 ... dV^iXtVtm, X k V k n k - 1) 

= (I[0(Xn i+1 (t)-XnM)I[®nM) (41) 
\ i=l i=l I 
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The conditions for forming the i th aggregate before t are the same as (|T3|), 
i.e. X^.(t) > X™i~ r (t), r = 1,2, ... ,rii — 1, where X^.(t) is the position of 
the subcluster of the first r particles in the i th cluster 

KM = Ui + ■ ■ ■ + n,_ x + ^tl) a + - £ v.t 



r s=l 



t 2 



+ jm(ni - r + n i+1 H h n k - m Ui ~^~2 ( 42 ) 

Here the velocities of the initial particles in this cluster have been simply 
labelled v±, v 2 , ■ ■ ■ , v ni and V ni = Z)™=i v s- Hence, as in (|2lD , the constraint is 

m— 1 



©«,(*) = J{e[x r ni {t)-xzc{t) 

i=i 

= n 1 4E^--K i + ^V^mr'] (43) 



i=l Vs=l n « 



We see from ([Tj]) that the factor 
fc-i fe-i 



n ec^c*) -x n M) = n « ( — - — ^ 

i=i »=i V n *+i ni * J 



in ([41]) depends only on the initial center of mass velocities V n Jrii of the 
clusters. Since the initial velocity distribution factorizes, we can perform the 
integration independently for each cluster, except for the variables V ni /rii 
that are coupled through ([44]) . Introducing for each cluster the change of 
variables ( p2|) one obtains as in (pT|)-(^) that the i th cluster contributes to 
the total integration on velocities in ( |4"1~D as 

J dvi . . .dv ni (f)x{vi) ■ ■■(j)x(v ni ) n [ v s ~ —V ni H rnrj . . . 



i=l \s=l 



= P n< ( r )_L- / dV ni exp • • • (45) 

Taking (|44|) and (^) into account in ([41]) and setting Ui = V rH /\^/n~i leads 
to the final result 



fJ-k (ni,n 2 , 
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x fdU x ... dU^{U x ) . . . <j>(U k ) n ( - -p= - "* t mr) (46) 

The first factor is the probability of formation of independent aggregates, 
whereas the second factor represents the correlations introduced between 
them by the gravitational forces. 

We now draw some important conclusions from the formula (|46D . We 
say that the i th aggregate is macroscopic if the number of its constituents 
Hi = T)N, < T] < 1, is a non vanishing fraction 77 of the total number of 
initial particles as N — > 00; its mass is then Mj = mrij = r\M tot . 

(i) Macroscopic mass 

According to the form of the mass density fl35|) at t = t* , it is clear that 
there can be only one macroscopic aggregate for t > t*. Let us recover this 
result using (f46|). The arguments of the ^-functions in (fi6|) are denoted by 

Ui+i Ui rii + n i+1 
w i>i+1 = — — — — mr 47 



The integral in fl46 ) is carried out on the domain T> of the variables Ui 
defined by u^j+i > 0. Suppose that the masses of the two clusters j and 
j + I become macroscopic, i.e. rij — > 00, — > 00 with mn^ = Mj > 0, 
mn,j + £ = Mj + £ > and rij remains finite for i 7^ j,j + £ Then as iV — > 00, 

w " +1 = V0rrit T < 48) 

»^ W =--^=-%*t (49) 

The relation ( |48l) together with the conditions iu^+i > for t > imply 
C/j + i > 0, . . . , Uj+£-i > 0, but (pE9|) implies also t/j+^_i < 0; thus the integra- 
tion domain X> shrinks to zero as iV — ► 00 and the corresponding probability 
vanishes. The argument is the same when more than two masses become 
macroscopic. 

We calculate now the probability to have one macroscopic mass, say rij = 
N — Z)i=y n ii holding the other aggregates rii, i ^ j finite, defined by 

)Ufc(ni, . . . , %_i, M tot , n j+1 , ...,n k ;t) = 
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lim n k (m, 



N ~ ^2 n h n j+i, ■■■,n k ;t) 

i¥=3 



(50) 



Notice that in the macroscopic limit, the mass M nj = M tot — mY^l^ni 
becomes infinitesimally close to the total mass: M n . — ► M tot . We make two 
observations on the probability P n (r), r > 0. If n is fixed 



lim P n {r) = P n (0) = - 



since this amounts to let m = M to t/N — > in 
fixed integer 



lim P N „ q (r) = P(t) 

N^oo 



(51) 

If n = N — q with q a 
(52) 



where 



P(r) 
with A(r) 



exp(-A(r)), 

OO 1 



2E 

n=l 



n J Mist T 



(53) 



is the probability of merging of the total number of particles into the single 
mass M tot . This is precisely the function determined in the proposition found 
in section 5 of ref. [4], because again m(N—q) = M tot —mq — > M tot as iV — > oo 
0. Finally, as iV — > oo the arguments u^j+i P7|) tend to 













M to4 




2A 




M tot 




2A 



T 



T 



(54) 



When (51]), (|52|) and ([54]) are taken into account in (|46|) (changing also f/j 
into — C/A one finds that the limit (I5DD is 



(55) 



H k {n X i nj_ x , M tot , rij+x, . . . , n k ; t) = 
Qj-xinx, 7ij_i; r)P(r)Q fc _ j (n fe , . . . , n j+1 ; r) 



2 In this proposition the quantity M tot /2A was set equal to one. 



16 



with 

Q i _i(n 1 ,...,nj-;r) = f dU x . . . dU^Ux) ■ ■ ■ <f>(Uj) 

ni ■ ■ ■ rij J 

xt (^-m... t (3=L.Jk.)„(3..^ T ) (56) 



The interpretation of ( |56|) is clear: the factor Qj-i (Qk-j) is the probability 
to find, left (right) of the macroscopic mass, j (k — j) aggregates made of a 
finite number of initial masses, that we call now fragments. We shall show 
in paragraph (ii) below that the probabilities ( |55| ) sum up to one. Therefore 
the only configurations that can occur after t* consist of a single macroscopic 
mass together with a dust of such fragments. One should notice that non 
macroscopic pieces of matter of the order mN v ', < u < 1, do not appear 
after t* . 

(ii) Statistics of fragments 

The probability to have a configuration of exactly k bodies after t* (i.e. 
one macroscopic mass and k — 1 fragments) is 

k oo 

M*) = X! H fJ>k(ni, ■ ■ ■ ,nj-i,Mtot,n j+ i, . . . ,n k ;t) (57) 

j=l n 1 ,...,n j - 1 ,n j+1 ,...,n k =l 

After the change of variables = (-^= ~ ^jf" t), one finds from (|56|) 

oo 

Qj-i( n i,---, n jl T ) (58) 



roo J 00 I / 

/ > d yi ...dy 3 ^—^f^ 

Jyi>-yj-i>yj>0 i=1 n=1 y/n \ 



where A(r) is the function defined in (|53f) . Hence from (|55|) , (|57|) , 
33!) one obtains the result 



and 



P(t) 



'A(t) 



fc-i 
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(Mr)) 



k-l 



exp(-A(r)), t>t* 



(k-l)\ 

As claimed above, the satisfy the normalization relation 



(59) 



^2vk(t) 



fc=i 

The distribution of the number k — 1 of fragments is Poissonian for all times 
t > t*. We find therefore that A(t) appearing in fl53|) has the interpretation 
of the mean number of fragments. A(r) tends to zero as a Gaussian when 
t —>■ oo and diverges as —2 In {^f- rj when r — ► (see eq. (52) in [4]). 

As a particular exemple, we write down from Q5"5D and (|5H) the probability 
of survival of a fragment of size n 



^ 2 {M toU n-t)= l -P{r) fl-erf 



lol 



Mtot 



2A 



r 



n(MtotT) 2 



if r 
if r 



(60) 



\J-KnMt 

where erf is the error function. 

The general position and velocity dependent distributions (40) for 
t > t* can also be written down more explicitely in the macroscopic limit. 
In particular, one has fi^X^Mtot; t) = 5(Xi - ^)5(V 1 )P(t). The macro- 
scopic mass is found at rest at the position M tot /2p without fluctuations, 
as expected. More generally, there are Gaussian small probabilities to find 
fragments far away from this point. 

So far we have given a full description of the state after t*: here the 
structure is simple since the weight of typical mass configurations is given by 



the set of distributions (45) with k finite. When t < t* the situation is more 
complex. Indeed, we have 



A' 



lim /if (ni,n 2 , 



n k ;t) = 0, k = 1,2, 



t < t* 



(61) 



since necessarily at least one of the rii tends to infinity and we know then 
from (ID that for r < 0, P n X T ) — V n ~~ * ( an the other factors in fldlf) 
are bounded by 1). Hence there remain always infinitely many aggregates 
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as iV — > oo, and the weight of typical configurations will be given by the 
distributions ^(ni, . . . , n&; t), k — > oo, involving infinitely many bodies. 
Computer simulations indicate that after a short transient time, typical con- 
figurations consist of approximately y/N aggregates, each of them having a 
mass of the order my/N. Thus we may conjecture that the distributions in 
this range, i.e. Y^^ n _ c • • • > n k\ t) with k ~ y/N, should have a 

non vanishing limit as iV — > oo. 

As far as the density (|20|) is concerned, we anticipate that it converges for 
t < t* to an absolutely continuous function, namely the uniform density (§) 
as in the static model. This would be consistent with the numerical observa- 
tion that after a short time during which most of the initial kinetic energy is 
dissipated by inelastic collisions, the subsequent evolution is dominated by 
the gravitational forces. The result of simulations is discussed in the next 
Section. Analytic proofs of these conjectures would complete the study of the 
dynamical phase transition that occurs at t* between a spatially extended 
and an aggregated phase of matter. 



4. NUMERICAL SIMULATIONS 



In this section we present the results of numerical simulations performed 
to determine the rate of formation of macroscopic masses for a system with 
a finite number N of particles evolving according to the model described 
above. In particular, we analyze the scaling to the continuum limit of the 
probability P^{t) for the formation of a macroscopic mass T]Nm, < r\ < 1, 
before time t. We also study the time evolution of the kinetic energy. 

Numerical simulations on this model are particularly simple, compared 
with their counterparts in higher dimensions, because the equations of mo- 
tion can be analytically integrated between succesive collisions. The simu- 
lation then reduces to keeping track of the particle masses, coordinates and 
velocities created in succesive collisions. 

In going to the continuum limit while keeping the density constant, dif- 
ferent scalings of the initial conditions are possible. Here, for reasons of 
numerical accuracy, we increase the number of particles N, while keeping 
the distance between them constant: a = 1, and also putting m — 1, so 
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that p = M/L = nm/na = 1. Further, we choose initial velocities with a 
Gaussian distribution of variance A = N/2. This procedure is equivalent to 



In Fig. 1 we show the mass formation probability P^{t) for iV = 1000 as 
a function of time, averaged over 10000 initial configurations, with 77 = 0.1, 
0.2,..., 0.9, 0.99, 1.0 from left to right. 

Note that P^if) clusters around the Heaviside function for 77 < 1, while 
the probability of total mass aggregation 77 = 1 follows the separate limiting 
curve given by (|53|). 

In order to study the scaling of Pj^it) curves as a function of the particle 
number N, we obtained the results for the initial values N = 10 • 2 r , r = 
0,1,..., 9, averaged over 1000 initial configurations. Fixing 77 = 0.5, we 
display in Fig. 2 the probability curves for increasing numbers of initial 
particles. They tend to the Heaviside function as iV — > 00. 

Having in view a quantitative study of this scaling we consider the time 
deviation \fl(N) — t*\ for increasing N, where t v ^(N) is the time at which 
the probability Pjy acquires the value /3j, i.e. P^(tp(N)) — j3, [3 = 
0.25,0.5,0.75,77 = 0.5. 

Fig. 3 reveals the power law 



The same behavior holds also for other values of 77. 

In order to get a deeper understanding of the dynamics of aggregation 
we analyze the evolution of the kinetic energy Ej^ n (t). If we started from a 
static initial configuration, we would simply find a parabolic law 



In Fig. 4 the ratio Ej? in (t) / E ki '° a (t*) has been plotted. This ratio approaches 
the normalized parabola as N — > 00. One can look at the local minimum of 
the curves as corresponding to the time scale of almost total dissipation of 
the kinetic energy due to initially numerous inelastic collisions. For N — > 00, 
the location of the minimum approaches the initial moment t = according 
to the power law iV -1 / 4 . Hence in the continuum limit the system gets 




1 




E, 



(t) = ^MU 2 
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instantaneously cooled down and the subsequent evolution is dominated by 
gravity. 

We also observe that at the time when the kinetic energy attains its 
minimum: 

(i) the average size of the formed masses scales as ~ m\fN 

(ii) the velocity distribution remains Gaussian with an effective standard 

deviation \ e ff ~ AiV -1 / 4 . 

It strongly suggests that before t* the density of mass in the continuum limit 
should also coincide with that of the static model (H). 

Finally, in order to study the sensitivity of our results to variations of 
the initial distributions, we have led the same simulations for interparticle 
distances following a Poissonian distribution and a Gaussian initial velocity 
distribution. We confirmed that all the results presented above remained 
identical, indicating a certain generality of the initial kinetic energy dissipa- 
tion process. 
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Fig. 1: Mass formation probability P£(t), N=1000, r) = 0.1, 0.2,..., 0.9, 
0.99, 1.0 from left to right. 

Fig. 2: Mass formation probability P^ 5 (t), iV = 10 • 2 r , r = 0, 1, . . . , 9, 
from left to right. 

Fig. 3: Log-log plot of \t}{N) - t*\, for (3 = 0.25,0.5,0.75 from up to 
down. 

Fig. 4: Normalized kinetic energy as a function of time for 
N = 10 • 2 r , r = 0, 1, . . . , 9. The parabolic limit curve represents the normal- 
ized kinetic energy of the continuous static initial configuration. 
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